Estimating von-Karman's constant from Homogeneous Turbulence 
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A celebrated universal aspect of wall-bounded turbulent flows is the von Karman log-law-of-the- 
wall, describing how the mean velocity in the streamwise direction depends on the distance from the 
wall. Although the log-law is known for more than 75 years, the von Karman constant governing 
the slope of the log-law was not determined theoretically. In this Letter we show that the von- 
Karman constant can be estimated from homogeneous turbulent data, i.e. without information 
from wall- bounded flows. 
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Introduction. The theoretical understanding of wall- 
bounded turbulent flows lags behind homogeneous turbu- 
lence, in which a number of universal constants and expo- 
nents can be estimated rather accurately on the basis of 
approximate arguments of considerable success [1] . One 
glaring such example of a lack of theoretical power con- 
cerns the apparently universal log law-of-the-wall which 
was discovered by von Karman in 1930 [2, 3]. The law 
pertains to the mean velocity profile V^(t) = ^Viyy) in 
wall bounded Newtonian turbulence (x, y and z are unit 
vectors in the stream-wise, wall-normal and span-wise 
directions respectively). In wall units the law is written 
as 

= — ln?;+ +B, for 30 < y+ < iicA . (1) 

Here the Reynolds number T^ba, the normalized distance 
from the wall and the normalized mean velocity 
V'^{y~^) (which is in the x direction with a dependence 
on y only) are defined by 
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Here p' be the fixed pressure gradients p' = —dp/dx, //q 
the kinematic viscosity. The law (1) is universal, inde- 
pendent of Tiex, the nature of the Newtonian fiuid and of 
the fiow geometry over a smooth surface, providing that 
Hex is large enough. 

It is one of the shortcomings of the theory of wall- 
bounded turbulence that the von Karman constant kk = 
0.44 ± 0.03 and the intercept i? 6.13 are only known 
from experiments and simulations [2, 4]. In this Letter 
we propose that kk can be estimated using universal con- 
stants that appear in homogenous turbulence. As such, 
we can draw on the relative power of homogeneous tur- 
bulence theory to improve our understanding of the char- 
acteristics of wall-bounded flows. We will not draw on 
any experimental information about wall-bounded flows. 

In constructing our argument we rely heavily on known 
facts, including recent results, concerning homogeneous 
isotropic turbulence and homogeneous anisotropic tur- 
bulence with a constant shear: dV{y)/dy — S. Due to 
Galilean invariance the statistics of the turbulent velocity 
field u{r, t) (from which the mean {u{r, t)) is subtracted) 



are independent of position. In other words, in such ho- 
mogeneous and anisotropic ensemble all statistical object 
computed at any point, like the density of the kinetic en- 
ergy K and the Reynolds stress W, 

X= (|M(r,t)|2)/2, W = -{u,ir,t)uyir,t)), (3) 

are space independent. On the other hand, two-point 
correlation functions, like the second order longitudinal 
and transverse structure functions, 

S2ir) ^ ((u||(r' + r,t)-«||(r',i))'), (4) 

S2{r) ^ ((«x(r' + r,t)-ru(r',t))'), (5) 

depend only on the vector separation r. In Eq. (4) and 
(5) M|| and uj_ are components of m, parallel and orthogo- 
nal to r. The physical reason for the homogeneity of the 
turbulent statistics is that the energy fiux e generated by 
the pressure head 
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is independent of the mean velocity itself (again, due to 
Galilean invariance) and is determined only by the space 
independent shear S and Reynolds stress W. 
I. Similarity of wall-bounded and constant-shear 
turbulence. We base our argument on the realization 
that wall-bounded turbulence in the log-law region and 
constant-shear turbulence are very similar. To be more 
precise, constant-shear homogeneous turbulence serves as 
a very good approximation to wall-bounded turbulence; 
various characteristics of turbulent statistics appear to 
coincide in the two flows within the available accuracy 
of physical experiments and numerical simulations. The 
basic reason for this similarity is precisely that the rate 
of the energy production (6) depends on the shear itself 
and not on its space derivatives, which obviously differ 
in these two flows. The cenral point of this Letter is that 
the similarity between wall-bounded and constant-shear 
turbulence allows one to estimate the von-Karman con- 
stant for wall-bounded turbulence using information from 
homogeneous constant-shear turbulence. 

The first result that we quote is long standing, stating 
a universal relation between W and K in a constant-shear 
fiow. 



W/K EE W+/K+ = 4 , CN « 0.53 



(7) 
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The the same value of cn, within the available accuracy 
(of about 5%), is measured also in the outer layer of 
channel flow [7]. This serves as additional support for 
the similarity between these two types of turbulence. 

Next we use the exact balance equation of mechanical 
momentum in a channel geometry, 

VoS{y) + W{y)=p'{L~y) . (8) 

In wall units (2) this equation reads 

S+{y+) + W+{y+) = l-y+/-RB^. (9) 

For 'Re\ ^ 1 the mean shear 6"+ in the log-law region (1) 
is governed by 

S+(t/+) = l/(«Kt/+), (10) 

For large values of y+ such that t/+ <C Vex Eq. (9) reduces 
to 

W+ = l, (11) 

meaning that in the log-law region the total momentum 
flux toward the wall is independent of the distance to the 
wall and is entirely accounted for by turbulent fluctua- 
tions (i.e. the Reynolds stress). Equations (7) and (11) 
show that in this region the kinetic energy is independent 
of the distance to the wall and has a universal value 

K+ = c^K (12) 

This constant value is shared by wall-bounded and 
constant-shear turbulence. 

II. Anisotropy in wall-bounded and constant- 
shear turbulent flows. The knowledge of the total 
kinetic energy density (12) is not sufficient for our pur- 
poses, we need to know the distribution of K between 
the different components of turbulent velocity, i.e. the 
values of 

K,^{ul)/2, Ky^{ul)/2, K,^{ul)/2. (13) 

The anisotropy of turbulent boundary layers, charac- 
terized by the dimensionless ratios Kj/K, plays an im- 
portant role in various phenomena and was a subject of 
experimental and theoretical interest for many decades, 
see, e.g. [2, 3]. Nevertheless, up to now the dispersion 
of results on this subject is too large. There is a widely 
spread opinion, based on old experiments, that the wall- 
normal turbulent fluctuations Ky are much smaller than 
the other ones. For example, in the classical textbook 
by A.S. Monin and A.M. Yaglom [2] it was reported that 
in a neutrally stratifled log-boundary layer = 0.54 ii", 
Ky = 0.06 if and if^ = 0.40 if. This appears to be 
in contradiction with simulation results at the largest 
Rca, 'Re\=590, which are available in Ref. [7]. These re- 
sults are reproduced in Fig. 1. Note the region of y+ 
values, about 100 < y+ < ^Vex, where the plots of 



Kj/K arc nearly horizontal, as expected in the log-law 
region. From these plots we can conclude that is this re- 
gion Ri 0.53 K which is close to the ratio 0.54, stated 
in [2]. In contrast, the simulational data for Ky/K and 
Kz/K are completely different. From Fig. 1 one gets 
Kj K. 0.22 K and K^, k. 0.27 K. Roughly speaking Ky is 
almost equal to K^- We propose that the difference be- 
tween Ky and Kz, observed in Fig. 1, is due to the effect 
of the spatial energy flux, that is expected to vanish in 
the asymptotic limit Tiex oo, but is still present at Tiex 
values available in direct numerical simulations [7]. In- 
deed, for both Vex shown in Fig. 1 Ky = K^ in the center 
of the channel, where the energy flux vanishes by symme- 
try. Clearly, there is no spatial energy flux also in a ho- 
mogeneous constant shear turbulent flow. Based on the 
similarity between the wall bounded and constant shear 
turbulent flows we expect the values of K^ to be the 
same for both flows in the limit Vex — * oo. The expec- 
tation Ky = Kz is confirmed by Large Eddy Simulation 
(LES) of a constant shear fiow [8]. As one sees in Fig. 2 
in this flow K^ w 0.46 /C, while KyP^Kz^ 0.27 K. For 
both flows Kz « 0.27 K, while in the stream- wise and 
and the wall-normal directions there is some difference 
in the ratios K^/K (0.53 vs. 0.46) and Ky/K (0.22 vs. 
0.27). We believe that these differences are again finite 
Tiex effects. This viewpoint is supported by a recent lab- 
oratory experiment by A. Agrawal, L. Djenidi and R.A. 
Antonia [9] in a vertical water channel with Re a = 1000 
which is reproduced in Fig. 3. Indeed the experimental 
values of K^/K, Ky/K and K^/K distribute similarly 
to channel simulations and constant-shear LES: 

K, = K/2, Ky = Kz = K/A. (14) 

The operational conclusion of this discussion is that 
there are experimental and simulational grounds to be- 
lieve that for both constant-shear and wall-bounded tur- 
bulent flows in the log- law region the turbulent kinetic en- 
ergy is distributed in a very simple manner: the stream- 
wise component contains a half of total energy, the rest 
is equally distributed between wall-normal and cross- 
stream components. We propose that Eq. (14) is asymp- 
totically exact. 

III. Structure functions in isotropic, constant- 
shear, and wall bounded flows. In isotropic homo- 
geneous turbulence the second order velocity structure 
functions, S2{r), [Eq. (4)] and S2{r) [Eq. (5)] are very 
well studied. They are invariant to the direction of the 
separation vector r and up to intermittency corrections 
they read 

S2{r) = C2(er)^= , &(r) = C2{er)^' , C2 = 2/3 . 

(15) 

The ratio of the dimensionless constants C2/C2 = 4/3 

follows from the incompressibly constraint, while the val- 
ues of C2, C2 are known from extensive experiments and 
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FIG. 1: DNS profiles of the relative kinetic energies in the 
streamvise direction, K^/K, lines denoted wall- 
normal, Ky/K, and spanwise, Kz/K, directions denoted as 

" " „„J li " l; T7„ cnn J^j- l; 
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FIG. 3: Experimental profiles of the relative kinetic energies 
in the vertical water channel with T(e\ = 1000 according to 
Ref. [9]. Solid lines serve only to guide the eye. Dashed lines 
show the energy distribution (14). 




FIG. 2: Relative components of the kinetic energies Ki/K for 
for constant-shear flow. Results of the LES [8]. Lines serve 
only to guide the eye. 

simulations [1, 2]: 

C2 = 4^2/3, C2«2.0. (16) 

For r larger than the outer scale of turbulence, the 
correlations between velocities in r-separated points van- 
ishes and the structure functions (15) saturate at their 
asymptotic values: 

S2{r)^2{u\), S2{r)^2{ul). (17) 

It is useful to define crossover scales I2 and £2 for the lon- 
gitudinal and transversal structure functions as follows: 

52(^2) = 2(uf) , S2{l2) = 2(tii) . (18) 



Clearly, in isotropic turbulence = (u^) and there- 
fore the scales £2 and £2 are related: 

?2 = {C2/C2f'^l2 = 3V3^2/8 « 0.65^2 . (19) 

The issue of the structure functions in anisotropic tur- 
bulent flows is considerably more involved [10, 11]. In 
addition to the isotropic contribution (15), the structure 
functions have anisotropic components, 5*2^-, belonging 
to irreducible representations of the SO (3) group with 

00 j 

S2{r)^J2 E (20) 

The various anisotropic sectors exhibit i dependent 
scaling exponents: 

5]^(r)(xr^fV,„,(r/r), (21) 

where the angular behavior of S'j^(r) is carried by the 
spherical harmonic Yj„i{r/r). The leading correction to 
the isotropic sector (15) is given the j = 2 contribution 

(2) A 

with Q « g. The mixture of contributions with dif- 
ferent different exponents, each with an amplitude that 
depends on the distance to the wall, may give the^ false 
impression that the scaling exponents of S2 and S2 are 
different; this impression disappears once the structure 
functions are projected on the various sectors of the sym- 
metry group, where their exponents are the same (ap- 
pearing universal) , but their amplitudes are of course 
non-universal, and see full details in [11]. 

Importantly, all the anisotropic contributions (21) can 
be eliminated by averaging over all the directions of r 
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because 

J Yjm {r/r) d{r/r) = , for j ^ . (22) 

After the elimination of the anisostropic sectors, the 
isotropic parts of 5*2 (r), ^2(^'), i-e. S'2,o(^'), 'S'2,o(^') ex- 
hibit the same scaling behavior as structure functions 
in isotropic turbulence, even in relatively low TZex chan- 
nel flows [12]. The situation for high TZe\ constant-shear 
and wall-bounded turbulence is even simpler. As demon- 
strated by Eq. (14), in these flows only the streamwise 
direction x is special, the partial kinetic energies in two 
other directions are equal: Kz = Ky. Therefore in the 
limit TZex oo the symmetry of these two flows can be 
considered as axisymmetric with the axis in the x direc- 
tion, and one can eliminate the anisotropic contributions 
by averaging only in the (x, z)-plane, parallel to the wall. 

This idea is supported by experiments in the atmo- 
spheric turbulent boundary layer [13]), in which the 
ground normal velocity component Uy was measured at 
point of fixed height y above the ground, separated by 
p which was parallel to the ground. Denote the result- 
ing structure function as S2{y,p)- After averaging over 
the azimuthal angle (j) in the (x,z)-plane, this function 
exhibits homogeneous scaling behavior (15) for all scales 
up to separations p in the (x, z)-plane which are close to 
the wall distance y: 

(S2{y,p))^ = C2{epf/\ p<y. (23) 

The constant C2 here is the same as in isotropic tur- 
bulence. Inspired by this experiment (and see Fig. 3 in 
Ref. [13]) we make the assumption that the largest of the 
two crossover scales (18), namely I2, is determined by 
the distance to the wall: 

i2=y. (24) 

In other words, the assumption is that the structure 
function (^S2{y,p)'^ changes sharply from a scaling law 

in p to its asymptotic constant value precisely at p = y. 
This assumption introduces an unknown factor of the 
order of unity to our arguments; we take this factor to 
be exactly unity. 

Relationship between C2, cn and kk- We have pre- 
sented all the ingredients necessary to estimate the von 
Karman constant. This constant will be related to the 
Kolmogorov constant C2 which appears in homogeneous 
isotropic flows and to the ratio cn which appears in ho- 
mogeneous constant shear flows. Using Eqs, (14), (18), 
(19) and (24) in Eq. (23) one gets: 

K = 4Ky = 2 {uD = C2{cl2f'^ = C2{eyf'^ . (25) 



In wall units K+ = C2{e+y+f'^. Taking e from Eq. (6) 
and if + from Eq. (12) we have 

c^'' = C2 {S+W+y+f^ . (26) 
Together with Eqs. (10) and (11) this leads to the rela- 
tionship: 

«K= (cnVC^)' . (27) 

Using the experimental values Cn ~ 0.53 and C2 ~ 2.0 
we get kn ~ 0.42 in excellent agreement with the known 
value of this constant, kn ~ 0.44 ± 0.03. It should be 
stressed that in fact we have used only one assump- 
tion (24) about the cross-over scale of the structure func- 
tion. All the other input is taken from homogeneous data 
without any wall in sight. We propose that the numerical 
agreement with known value of kk reflects the quality of 
the input values of C2 and cn- 
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